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In this paper, we derive a general mixed (bright-dark) multi-soliton solution to a one-dimensional 
multicomponent Yajima-Oikawa (YO) system, i.e., the {M-\- l)-component YO system comprised 
of M-component short waves (SWs) and one-component long wave (LW) for all possible combina¬ 
tions of nonlinearity coefficients including positive, negative and mixed types. With the help of the 
KP-hierarchy reduction method, we firstly construct two types of general mixed A^-soliton solu¬ 
tion (two-bright-one-dark soliton and one-bright-two-dark one for SW components) to the (3-1-1)- 
component YO system in detail. Then by extending the corresponding analysis to the {M-\- 1)- 
component YO system, a general mixed V-soliton solution in Gram determinant form is obtained. 
The expression of the mixed soliton solution also contains the general all bright and all dark N- 
soliton solution as special cases. Besides, the dynamical analysis shows that the inelastic collision 
can only take place among SW components when at least two SW components have bright solitons 
in mixed type soliton solution. Whereas, the dark solitons in SW components and the bright soliton 
in LW component always undergo usual elastic collision. 
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1. Introduction 

The investigation of multicomponent nonlinear systems has received much attention in recent 
years. Of particular concern is the multicomponent generalization of the nonlinear Schrodinger 
(NLS) equation, namely, the vector NLS equation.Particularly, it has been shown that mul¬ 
ticomponent bright solitons undergo shape changing collisions with intensity redistribution.^”^ ^ 
This interesting collision feature has been widely studied in real physical systems such as nonlin¬ 
ear optics and Bose-Einstein condensates. 

The long-wave-short-wave resonance interaction (LSRI) is a fascinating physical process, in 
which a resonant interaction takes place between a weakly dispersive long-wave (LW) and a short¬ 
wave (SW) packet when the phase velocity of the former exactly or almost matches the group 
velocity of the latter. The theoretical investigation of this LSRI was first done by Zakharov^^ on 
Langmuir waves in plasma. In the case of long wave propagating in one direction, the general 
Zakharov system was reduced to the one-dimensional (ID) Yajima-Oikawa (YO) equation.This 
model equation also appears in diverse areas that include hydrodynamics,^^ nonlinear optics,^®’^^ 
biophysics^^ etc. The ID YO system is integrable by the inverse scattering transform method^^ 
and admits both bright and dark soliton solutions.The rogue wave solutions to the ID YO 
system have recently been derived by using the Hirota’s bilinear method^^ and Darboux transfor- 

mation.26’27 

In the present paper, we consider a general multicomponent ID YO system^^ 

i5f)-5S+L5W=0, f=l,2,---,M, 

L, = 2LtiO^|sW|2, o^ = ±l, 

which describes the dynamics of nonlinear resonant interaction between a LW (L) and multiple 
(say M) SWs (S^^^) in one-dimensional case. Hereafter, we refer to the above multicomponent 
system as the (M-i-l)-component YO system. It was pointed out that the (2-i-I)-component YO 
system can be deduced from a set of three coupled NLS equations governing the propagation of 
three optical fields in a triple mode optical fiber by applying the asymptotic reduction procedure.^^ 
Eqs.(l) has been derived to describe the interaction between a quasi-resonance circularly polar¬ 
ized optical pulse and a long-wave electromagnetic spike.In the context of the many-component 
magnon-phonon system, such multicomponent YO-type system has also been proposed and its 
corresponding Hamiltonian formalism was studied.^^ Also, the authors in ref.^^ have carried out 
Painleve analysis for Eqs.(l) and obtained the general bright N-soliton solution in the Gram de¬ 
terminant form. Later on, they constructed an extensive set of exact periodic solutions in terms of 
Lame polynomials of order one and two.^^ The rogue waves of the (2-1-1 )-component YO system 
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with Oi = 02 = 1 have been reported in Ref..^^ 

It is worth mentioning that the KP-hierarchy reduction method to derive soliton equations as 
well as their solutions was developed by the Kyoto school.A number of soliton equations such 
as the NLS equation, the modified KdV equation, and the Davey-Stewartson equation can be re¬ 
duced from the KP-hierarchy. Indeed, the multicomponent YO system (1) with same nonlinearity 
coefficients (all O/ = l,or O; = — 1) was derived from the KP-hierarchy in.^^'^^ In particular, 
general (pseudo-) reductions of the two-dimensional Toda lattice hierarchy to constrained KP sys¬ 
tems with dark soliton solutions were introduced in^* and reductions to constrained KP systems 
with bright soliton solutions from the multi-component KP hierarchy were introduced in.^^ For 
the (l-i-l)-component YO system (when M = 1 in Eqs.(l)), the detailed interpretation of how to 
obtain the Wronski-type bright and dark soliton solutions by using the reduction technique was 
presented in Ref..^® By applying this method, general dark-dark soliton solution was derived in 
a coupled focusing-defocusing NLS equation.^^ Recently, one of the authors^^ has constructed 
general bright-dark A-soliton solution to the vector NLS equation of all possible combinations of 
nonlinearities including all focusing, all-defocusing and mixed types. In a previous study,^^ we 
have obtained both the Gram-type and Wronski-type determinant solutions of N dark-soliton for 
Lqs. (1). 

The goal of this paper is to construct general bright-dark multi-soliton solution to the multicom¬ 
ponent YO system (1). The rest of the paper is arranged as follows. In section 2, we deduce two 
types of general mixed soliton solution, which include two-bright-one-dark and one-bright-two- 
dark soliton for SW components, to the (3-1-1 )-component YO system by using the KP-hierarchy 
reduction technique. In section 3, general bright-dark soliton solution consisting of m bright soli- 
tons and M — m dark solitons to the multicomponent YO system (1) is obtained by generalizing 
our analysis. We summarize the paper in section 4. 

2. General mixed-type soliton solution to the one-dimensional (3+l)-component YO system 

Prior to presenting the general mixed-type soliton solution for Eq.(l), we first consider the ID 
(3-i-l)-component YO system 

= 0 , 

i5;^^-5i?+L5(2)=0, 

i5j^^-5i?+LS(3) =0, 

Lt = 2(oi p + 0215(2) |2 + 0315(3) |2)^, 

where o^ = ±1 for f = 1,2,3. Lor three short-wave components, there exist two kinds of mixed- 
type vector solitons (two-bright-one-dark and one-bright-two-dark). Hence, in the subsequent two 
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subsections, we will construct these two types of soliton solutions, respectively. 


2 .1 Two-bright-one-dark soliton for the SW components 

In this case, we assume that the SW components 5*^^^ and are of bright type and the SW 
component is of dark type. By introducing the dependent variable transformations 

= = 5(3)_p^^gi(aix+a2f)^ L=-2(l0g/)«, (3) 

equations.(2) are converted into the following bilinear equations: 


[iZ),-D2]g«./ = 0, k=l,2, 

[i(A-2aiD^)-D2]/j(i)./ = 0, (4) 

(AA - 2o3p2) / ■ / = -2 LLi - 2o3plhWh(^>, 

where and h^^^ are complex-valued functions, / is a real-valued function, ai and pi (pi > 

0) are arbitrary real constants and * denotes the complex conjugation hereafter. The Hirota’s bilin¬ 
ear differential operator is defined by^^ 

S-a?) (a--3?) 

In what follow, we will show how the (3-1-1 )-component YO system and its mixed type mul- 
tisoltion solution can be obtained from three-component KP hierarchy by suitable reductions. To 
this end, we firstly present tau functions of the Gram determinant form in three-component KP 
hierarchy 
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where7 is an x identity matrix, A and B are N xN matrices whose entries are 
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and 0 is a N-component zero-row vector, 4>, 'P, T, 'P and T are /7-component row vectors given 
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$ = • • • ,e^^), 


vp = (e^i,e’i 2 ,-- 

•,e'i^), T=(e^i,e^^-- 

• ,e^^ 

vp = (efli,e^2 

• ,e^^), T = (e^ke^2, • • 

• ,e^^ 


with 

t" ^(1) 2 t-F ^(1) - -2 F 

Q = - Xl{ + PiXi + PiX2 + g-O, W --^_i + -PjX2 + ^j0, 

Pi-Cl Pj+Cl J J 

T]i = qiy[^^+T]io, ^j = qjy[^^+'^jo, Xi = riyf^+Xio, %; = 07f^ + %;o- 

Here PnPj,qi,qj,ri,rjXiO,^jO,^iO,^jO,XiO,XjO and ci are complex parameters. Based on the Sato 
theory of the KP hierarchy,one can find that the following bilinear equations are satisfied by the 
above tau functions 

(f)x2 -■D^i)^Lo(^l) ■T^O,o(^l) = 0, 

(^X2 -■Dxi)^0,l(^l) ■T^o,o(^i) = 0, 

{Dx 2 -DI -2ciDxi)xofl{ki + 1) ■Xo,o(fci) = 0, 

^xi^^a)Xo,o(^i) ■Xo,o(^i) = -2xi,o(fci)x_ip(fci), 

^xi^y( 2 )Xo,o(^i) ■Xo,o(^i) = -2xo,i(fci)xo-l(fcl), 

{DxiD ( 1 ) -2)xo,o(fci) ■Xo,o(^i) = -2xo,o(fci + l)xo,o(^i - !)• 

A_i 

Here we omit the proof of the above bilinear equations, which can be given by using the Grammian 
technique.^^’^^ Instead we carry out the reduction process to obtain bilinear form (4) of the (3+1)- 
component YO system. By setting xi, yf''^ are real, X 2 , ci are pure imaginary and by 

letting p* = pj,q*=qj = r* = Pj, ^*q = Ijq, ri*^ = fj jq and X*q = Xjo, one can check that 

aijih) =a*i{ki), bij = b*i. 


Furthermore, we can define 

/ = Xo,o(0), =Xi,o(0), g(2) =Xo,o(l), 

^(1)* = -x_po(0), = -xo,-i(0), = xo,o(-l), 


and thus the bilinear equations (8) become 

(D,,-Dj)g«./ = 0, 

(D,,-Dj-2ciD;,Jfi(i)-/ = 0, 

Dx,D,^f-f = 2g(^)g(^>, k=l,2, 
>1 

(D,,D(i)-2)/-/=-2fi(i)fi(i)*. 

^-1 


( 9 ) 
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Next, we conduct the dimension reduction by rewriting / as 

A' / 

/= , , 

-/ B' 

via performing row and column operations, where A' and B' are N x N matrices defined as 

a'. .(kA- __ b'- ■ - __eX,+x*H*Hy 

+ + r, + rf 


+ — ^iy[ * H ——Al +p*x\ — pyx 2 + ^*o+y]iQ, 

Pi + Cl 

T\j + + PFi + pj^2 + ^jo + ri*o, 

Xi + ^1 = fiy[ ^ H-TT--^“1 P*i^'^ ~ P^^-’C2+^*o + XiO, 

Pi “C 

%; + = r*jyf^ + + PjXi + p]x2 + ^jo + X)q, 

Therefore, under the reduction conditions 

*2 2 i^3p? 

Pi =ic>i^,-——, Pj = -\( 5 iqj + -, 
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1^3 Pi 1 ’ 
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the following relation holds 


and thus one can get 


^X2b'ij = (-ioid (1) -i023 (2) +i 03 pjd,c^^)bA 


fx2 = -i<5i/(1) - i<52/(2) +i03pi/x_j, 
21 21 


and its derivative with respect to jci. 


fxiX2 = “^^2/ (2) +i03pi/;i:j;c_l. 

•^1/1 -^1/1 


( 17 ) 
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On the other hand, the last three bilinear equations in (9) expanded as 


fmf-hfm =«<V. 

■^ 1/1 71 -^ 1/1 71 


and 


(18) 


give 




-ifxix2f + ^fx2fxi-03plf = 


(19) 


( 20 ) 


by referring to Eqs. (16)-(17). 

Applying the variable transformations 


Xl —X, X2 — —It, 


( 21 ) 


i.e., 


dx — dx ,, dt — idj;2, 


and taking ci = itti, the first three equations in (9) become the first three bilinear equation in (4) 
and Eq. (20) is nothing but the last bilinear equation in (4). 

Eastly, let e’^'O = c^*, ~ ~ (*' ~ I 52 , ,N) and define Ci = 

— (cj^\c 4 ’ ■ ■ ■ ’ 4 ^) ^2 = ~(c[^\c 4 ’ ■ ■ ■ arrive at the general mixed solition 

(two-bright-one-dark soliton for SW components) solution to the ID ( 3 - 1-1 )-component YO system, 

A 7 4)^ 
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where A, a(^( and 5 are A x A matrices whose entries are 
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and (|) and Q are A-component row vectors given by 
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with 0, = piX — ipjt + Qio and pi, 0,o and c] (/ = 1,2, • • • ,A;k = 1,2) are arbitrary complex 
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parameters. 


2.2 One-soliton solution 

To obtain the single soliton solution, we take = 1 in the formula (22). The Gram determinants 
read 


(23) 
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where an = ^e®>+®saS;^ = ^ (-ff+S) = 

for k= 1,2. Then one ean write the above tau funetions as 

/= i^QQi+&*i+2R{i,n^ 

^(1) =c5^^e®', 


•\-2 W* W 


iPl+PlM-Pl + H^rrSk- 
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(24) 

(25) 

(26) 


with 




{k)*(k) 


^Vk=iOkcr c 


ip\+piWi -P1+ (p|+i«oTJi-i«i) ] 


c — 


P\ + lai 


It is necessary to note that this mixed soliton solution is nonsingular only when 

The above tau functions lead to the one-soliton solution as follows 

ik) 

5^ = ^e^^(^’i*)e'®''sech[0i^ + 7?(l,r)], k= 1,2, 

5(3) _ p^gi(an+a2,)^^ - (1 -e2i^')tanh[ei/j -f 7 ?(l, 1*)]}, 

L=-2p2^sech2[0i^+7?(l,r)], 


> 0 . 


(27) 

(28) 
(29) 


where 0i = 01 /; + i0i/, the suffixes R and 1 denote the real and imaginary parts, respectively. The 
quantities =1,2 represent the amplitude of the bright soliton in the SW components 

and the real quantity 2pj^ denotes the amplitude of soliton in the LW component —L. For 
the dark soliton in the SW component approaches |pi| as x —)■ ± 0 °, and the intensity 

is |pi|cos(|)i. As an example, we illustrate the mixed one-soliton at time t = 0 in Fig. 1 for the 
nonlinearity coefficients ( 01 , 02 , 03 ) = (1,-1, !)• The parameters are chosen as pi = 1, tti = 2, 
Pi = I -l-i, 010 = 0, = 1 and (a) = 3 -|-i; (b) = 1 -|-i. One can observe that when 

the parameters take different values, the intensities of bright solitons for SW components 
and change, but the depths of dark soliton for SW component and of the soliton for LW 
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component L remain the same. 



X I 


Fig. 1. (Color online) Mixed one-soliton solution (two-bright-one-dark soliton for SW components) in (3 h- 1)- 
component YO system. 


2.3 Two-soliton solution 
By taking // = 2 in (2: 

t 

( 

f = 


where atj = < 

k= 1,2. Then the explieit 

/ = i+£(i,r)( 

+£(1,1*,2,: 


I, we have the following Gram 


bi 

a \2 1 

0 


111 

022 0 

1 


-1 

0 bn 

b\2 


0 

— 1 b2i 

b22 


an 

ai2 

1 

0 

an 

022 

0 

1 

-1 

0 

bn 

bi2 

0 

-1 

bii 

bii 

0 

0 - 

Ak) 

‘'1 

M 

k'l 


1 { 

Pi-iai 

\ gOi+e 

^ij ■ 

- Pi+P* \ 

p}+iai 

J® 


form of the above tau funetions 


. 01+91 +£( l , 2 *) e ®'+®2 +£( 2 , 1 *) 
7*^e0i+02+0i+02 

e02_pi7(*)(1^2,r)e®‘+®2+91_^^(^ 


ill 

«12 

1 

0 

^21 

«22 

0 

1 

-1 

0 

bn 

bi2 

0 

-1 

bii 

bii 


{Pl+Pj)[P*i -Pj+ (p*+ia;)] 

.n be expressed as 

;92+91_p£(2,2*)e®2+02 

(32) 

'(l,2,2*)e®i+®2+0^, (33) 
































J. Phys. Soc. Jpn. 


DRAFT 


+G(l,l*,2,2*)e®'+®2+6I+62, 


(34) 


where 
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Fig. 2. (Color online) Mixed two-soliton solution (two-bright-one-dark soliton for SW components) including in¬ 
elastic collisions for the SW components 5^^ and in the (3H-l)-component YO system . 






Fig. 3. (Color online) Mixed two-soliton solution (two-bright-one-dark soliton for SW components) including elas¬ 
tic collisions for the SW components 5^*^ and 5*^^^ in the (3H-l)-component YO system . 


The above solution eontains both singular and nonsingular solutions. To assure the nonsingular 
solution, the denominator / needs to be real and nonzero. The expression for / ean be rewritten as 

f = 2 e ®'*+2 cosh( 0 ii;- 02 i;H- - -)+e^'"eos( 0 i/- 02 / + //i/) 
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+e 2^cosh(0i^-e2i; + ^)], (35) 

where 

e^i=£(l,r), =£( 1 , 2 *). 

Combining the requirement for the existenee of one-soliton solution, we ean eonelude that the eon- 

^3 TJ 

dition£(z, r) >0,z = 1,2 is a necessary eondition and e 2 -\-q~ > is a sufficient condition 

to guarantee a nonsingular two-soliton solution. For the interaction properties of these solitons, one 
can carry out the asymptotic analysis as and deduce that the bright solitons appearing 

in SW components S^'\i =1,2) only undertake elastic collisions under some special parameters, 
while the dark solitons in the SW component and the bright soliton in the LW component —L 
always exhibit elastic collisions. More precisely, bright solitons in SW components =1,2) 

undergo elastic collisions if the parameters satisfy the condition ^ = -[jy- Otherwise, they under- 

^2 ''2 

take inelastic collisions (shape changing). For illustrative purpose, the interactions of two solitons 
are displayed in Figs. 2 and 3 for the (3-1-1 )-component YO system with the nonlinearity coefficients 
( 01 , 02 , 03 ) = (1, — 1,1). The parameters in Fig.2 are chosen as pi = l,ai = l,pi = | — |i,P 2 = 
l-|-^i,c[^^ = l-|-i,C 2 ^^ = 2 ,c[^^ = ^i,C 2 ^^ = 1 -|- 2 i and 010 = 020 = O 7 which result in inelastic colli¬ 
sions for the bright solitons in the SW components 5*^'^ (z = 1,2). Meanwhile, we show an example 
of elastic collision between bright solitons in two SW components in Fig.3 with the parameters 
pi = l,ai = l,pi = | - |i,p 2 = l + y^i,ci^^ = 2 , 4 ^^ = l,cf^ = 1 , 4 ^^ = I and 010 = 020 = 0. In 
Figs. 2 and 3, the SW components 5*^^^ and 5*^^^ display two bright solitons collisions in (a) and (b), 
the SW component shows two dark solitons collision in (c), and (d) represents two bright soli¬ 
tons collision of the LW component —L. As shown in (a) and (b) for the SW components collision 
of and respectively, one is elastic collision (2), the other is inelastic collision (3). 



Fig. 4. (Color online) Mixed two-soliton bound state (two-bright-one-dark soliton for S W components) in the (3+1)- 
component YO system. 


In addition, soliton bound states are one class of special multisoliton solutions, in which mul¬ 
tiple solitons move with the same velocity. By assuming pi = pm + ipu, one can obtain the mixed 
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two-soliton bound state from (32)-(34) with the restriction pu = p 2 i. To demonstrate such a bound 
state, we choose the parameters as pi = l,ai = l,p\ = \-\-\^,P 2 = 1 + = 2 + i,C 2 ^^ = 

2 + 4i, = ^ + i, ^ + 2i, 010 = 02o = 0 under the same nonlinearity coefficients as in Figs. 

2 and show the contour plots of all components in Fig. 4. Here, Fig. 4 (a) and (b) exhibit two- 
bright-soliton bound states for the SW components and with the different amplitudes. 
Fig. 4 (c) displays two-dark-soliton bound state for the SW component 5*^^^ and two-bright-soliton 
bound state of the LW component —L is shown in Fig. 4 (d). 


2.4 One-bright-two-dark soliton for the SW components 

In this subsection, we assume the SW component 5*^^^ is of bright type and the SW components 
and are of dark type. The dependent variable transformations 

5(1) = 5(2) =pi^ei(“>^+“^), 5(3) = p2^ei(“2^+“20, -l{\ogf),^, (36) 

convert the (3-1-1 )-component YO Eqs.(2) into the following bilinear equations: 

[iT>,-D2]g(i)./ = 0, 

[i(A-2a;D;,)-D2]/i(0./ = o, / = i,2, (37) 

(AA - 212^1 0;+ip2) / . / = -20ig(l)g(l)* - 212^1 O/+ip2/i(0/i(0*, 

where g(i),/z(i) and /z(2) are complex-valued functions, / is a real-valued function, a/ and p/(p; > 
0; / = 1,2) are arbitrary real constants. 

Here the bilinear form of the (3-1-1 )-component YO system (37) is viewed as a reduction of 
two-component KP hierarchy. To this end, we start with the tau functions expressed in Gram de¬ 
terminants as follows 


xo(^i, A) = 


A 7 
-7 B 


(38) 


(39) 


where d>, 'P, d>, 'P are /7-component row vectors defined previously, A and 5 are /7 x /7 matrices 
whose entries are 
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Tl.= 


T-^ 

Pj + Ci Pj+C2 

( 1 ) , - - ( 1 ) 

+11/0, ^j = (ijy\ 


+ PjXi -p']x2+^j0, 

+ TijO- 


Here and C 2 are complex parameters. Based on the Sato theory for 

the KP hierarchy,the above tan functions satisfy the bilinear equations 

{Dx2-Dl^)xi{ki,k2)-Xo{ki,k2) = 0 , 

{Dx 2-DI^ -2ciDx^)XQ{k\ + l,k2)-XQ{ki,k2) = 0, 

{Dx2-DI^ -2c2£>xi)Xo(fcl,^2+1) ■Xo(^1,^2) =0, 

Dx^D (i)XQ{k\,k2) ■XQ{k\,k2) = —2xi{ki,k2)x-\{ki,k2): 

y\ 

{DxiD (\) —2)Xo(fcl,fc2) ■Xo(fcl,^ 2 ) = —2Xo(fcl + l,fc 2 )Xo(^l — 1 ,^ 2 ), 

X _1 

(DxiD (2) —2)xo{ki,k2) ■Xo{ki,k2) = —2xo(fci,fc2 + l)'To(^i5^2 — 1), 

't_l 

Next we perform the reduction process to obtain the bilinear equations (37). We first consider the 
complex conjugate reduction by setting xi, to be real, X 2 , ci C 2 to be pure imaginary 

and by letting p*- = pj, q*j = qj, ^*q = ^jo, and ri*^ = fjjo- Then, it is easy to see that 

(^ij{ki,k2) = ^jii ki, ^ 2)5 kij = bjj. 


Thus, it then follows 

/ = Xo(0,0), g(i)=xi(0,0), fi(i)=Xo(l,0), 

^(i)* = -x_i(0,0), Mi)*=Xo(-l,0), fi(2)*=Xo(0,-l), 


the bilinear equations (40) become 

(D„-D?,)gl‘)./ = 0 , 


(D 


'Xl, 

X2 ^Xl 


2ciDx,)h(^^-f = 0, 


Dx,Dwf-f = 2 g(^)gW\ 

•>^1 

(Dx.D^m -2)/-/ = -2fi«fiW*, /=1,2. 

•^-1 

Similar to the two-bright-one-dark soliton case, we can show that if qi satisfies 


,,'2 = ir, „ _ ‘”2Pi _ ig3P2 

p* + Ci p; + C2’ 

J J Pj-Ci Pj-C2 


i.e.. 


101 


+ {p*+Pj)[p*-Pj + 

(v*+ci](n 


+ 


1^3 P 2 


(P*+Cl)(/3j—Cl) ^ (p*+C2)(pj-C2 


(41) 


(42) 

(43) 

(44) 
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one can get the following relation 

^2 = -ioi/(1) + i(52p?/(1) + iaspi/(2), (45) 

/I -^-1 -^-1 

and its derivative with respect to xi, 

fx,x2 = + i(53pi/^^^(2). (46) 

On the other hand, the last three bilinear equations in (41) expanded as 

/u)/-A,/i.) (47) 

■^1/1 /I 

and 

/ Mf-f.Jm -f = -A"'*'"*, fmf-fnfm-f = -A'"**®*. (48) 

AJA_A_j A]^A_j A_j 

give 

-i/xiX2/ + i/x2/xi - (02pf+ 03pi)/ = -03pi/z(^(/z^^^*, (49) 


by referring to Eqs. (45)-(46). Eq.(49) is exactly the last bilinear equation in (37) through the 
variable transformation (21) and ci = iai,C 2 = i0C2. Under the same transformation, the first three 
bilinear equations in (41) become the first three in (37). 

In summary, we complete the reductions from the bilinear equations in (40) to the ones in 
(37). Thus, we are able to construct general mixed soliton (one-bright-two-dark soliton for SW 
components) solution to the the ID (3-1-1 )-component YO system, 

(n ^ ^ 

^ I m I (u 

f= -I B 0^ , (50) 

-IB -IB 

0 Cl 0 

where A, a(^) and B are N xN matrices whose entries are 



(0 

a - ■ = 


bij = 


Pi + P*, 


Pi-^^i \ g9,+e* 
p) + iaij 

( 1 )* ( 1 ) 

lOic; ’ c)' 

^ J 


(pI + Pj) [p 1 - Pj + (p*+iai)(/7‘2-i(Xi) + (P*+ia2)(pj-mJ 

and (|) and Ci are /7-component row vectors given by 


(^ = (e®^e®^•••,e«^), Q =-(ci^^ 


(1) di) 




with 0,- = piX — ipjt -f 0/0 and pi, 0/o and c- = e’l'° {i= 1,2, • • • ,/7; / = 1,2) are arbitrary complex 
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parameters. 

Remark 3.1 In the last subseetion, we have eonstrueted the soliton solution in which the SW 
components 5*^^^ and are of bright type and the SW component 5*-^^ is of dark type in the (3+1)- 
component YO system. It is noted that we start from the (2+l)-component KP hierarchy with one 
copy of shifted singular point (c). In contrast, for the soliton solution (the SW components 5*^^^ is 
of bright type and the SW component and are of dark type) obtained in this subsection, we 
begin with the (l+l)-component KP hierarchy with two copies of shifted singular points (ci and 
C 2 ). Thus the number of the components in KP hierarchy matches the numbers of the components 
possessing the bright soliton while the number of the copies of shifted singular points coincides 
with the numbers of the components possessing the dark soliton. This fact can be refered to the 
construction in Ref..^^ 


2.5 One-soliton solution 

By taking = 1 in the formula (50), we get the Gram determinants 


/ = 


<311 1 


(0 

1 


ail 

1 

e®' 

, = 

a\l 


-1 

bn 

0 

-1 bn 


-1 

bn 


fi') 







0 

_ ^ 

0 


(51) 


where an = 


pi+p\ 




pi-iai \ pOi+e^ 


pi+p*i \ pl+i«/y 


1 and bn = 


iGiC 


( 1 )* ( 1 ) 

1_^ 


(p*+pi)Ai 


Ll 


iPf+lPF 


1 {p*,+iai){pi-iai) 


for I = 1,2. These tau functions can be rewritten as 

= /=1,2, 


-, All = [P*i-Pi + 

(52) 

(53) 


with 

( 1 )* ( 1 ) 

„ 2 i(|), _ Pi 

(p*+pi)2Air pi+iar 

2 

where An = [p\ — p\ +LLi (p*+ia;|(’pi-la;) ] • mixed soliton solution is nonsingular 

only when > 0. 

In this case, the one-soliton solution has the following form 

Ai) 

5 ( 1 ) = ^e^^(^’i*)e'®''sech[0i^ + R(l,r)], (54) 

= p,e‘(“'^+“''){l+e2‘^'-(l-e2i^')tanh[0i^-fR(l,r)]}, /= 1,2, (55) 

L = -2p2^sech2[0i^+R(l,r)], 


(56) 
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where 0i = diR + idu, the suffixes R and I denote the real and imaginary parts, respeetively. The 


quantity 


^( 1 ) 

represents the amplitude of the bright soliton in the SW eomponents 5*^^^ and 
the real quantity denotes the amplitude of the soliton in the LW eomponent —L. For the dark 
soliton in the SW eomponents and | approaehes |p/1 as x —>■ ±°o, and the intensity is 

IP/1 eos (|)/ for / = 1,2. When tti and a 2 take different values, there are two oases: (i) ai = a 2 . In this 
case, (|)i = (|)2 means dark solitons in SW components and are proportional to each other. 
Thus this situation is viewed as degenerate case for dark solitons. (ii) ai a 2 . The condition (|)i 
(|)2 implies that dark solitons in SW components and have different degrees of darkness at 
the center. In this situation, the SW components 5*^^^ and are not proportional to each other. We 
illustrate such degenerate and non-degenerate cases for the choice of the nonlinearity coefficients 
( 01 , 02 , 03 ) = (1,1, — 1) in Fig. 5. The parameters are chosen as pi = 1, pi = ^ — i, Qio = 0, = 1 

and (a) tti = 02 = —1, P 2 = \\ (b) ai = —|,a 2 = —P 2 = 1 at time t = 0. 




Fig. 5. (Color online) Mixed one-soliton solution (one-bright-two-dark soliton for SW components) in (3+1)- 
component YO system. 


2.6 Two-soliton solution 

When N = 2 m (50), the Gram determinants are of the form 


an 

an 

1 

0 


(0 

a\l 

“12 

1 

0 

an 

ail 

0 

1 


(i) 

a\{ 

(0 

«22 

0 

1 

-1 

0 

^11 

b\2 


-1 

0 

^11 

^12 

0 

-1 

b2i 

b22 


0 

-1 

^21 

bn 


( 57 ) 
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where atj = 



an 

an 

1 

0 

e®' 


an 

«22 

0 

1 

0^2 


-1 

0 

hit 

^12 

0 


0 

-1 

bn 

bn 

0 


0 

0 

- 

‘'1 

‘'2 

0 

1 e0,+0* JO _ 

1 

( Pi-'m 

4 e0,+e 

2 and 

Pi+pf - 

Pi+P* 

V P*j+'iai 



(58) 


( 1 )* (1 
iGiC) ’ c) 

* J 


^^=1 (p*+ia;Xp'-ia/) ] ^ 2. Then we ean express the tau funetions for two-soliton solution as 

/ = l+£(l,r)e®'+®^+£(l,2*)e®i+®2+£(2,r)e®2+®(+£(2,2*)e®2+®2 


+£(l,r,2,2*)e®i+®2+01+05^ 

^(1) = c[^)e®'+4^^e®2+F(l,2,r)e®‘+®2+0I+^('1^2,2*)e®'+®2+05^ 

= 1 + g(') (1, r )e®' + g(') (1,2*)e®' +®2 + (2,1*)e®2+®( + G^'^ (2, 2 *)e® 2+02 

+g('4i, l*,2,2*)e®'+®2+0I+05^ 

where 

G«(hr)- 


(59) 

(60) 

(61) 


E{iJl = 


( 1 )* ( 1 ) 
lOic; ’ c) ’ 

^ J 


igf+ipf 


{Pi+P*j)^[p--Pj + i:i=i (p*+ia,)(p,^icx,) J 
,r /2n 1* 

£(l,r,2,2*) = |pi-p2| 


Pj+m 


E{i,f 


|2 

[ £(l,r)£(2,2*) 

£(1,2*)£(2, r) ] 


JP\+P\){P2+P\) 

{Pi+P\){P2+pl)^ 


G<»(i.r.2.r)= <;;;;;g|g;;g £(i.r.2.r). 

Same as the previous subseetion, nonsingular solution requires the denominator / to be real 
and nonzero. For this purpose, we rewrite / as 

r, . a ^1+^2 — ^2? u 

/ = 2e®'*+®2fi[e^— cosh(0ip-02pH-r-)+e^'*eos(0i/-02/ + Ffi/) 


+e 2 ^eosh( 0 ip- 02 p + ^)], 


(62) 


where e^i =£(1,1*), = E{2,2*), ,2,2*), e^i«+‘-^i2 = £(i^2*). Thus, one 

ean eonelude that E{i,i*) > 0,i = 1,2 is a neeessary eondition to obtain a regular solution and 


^1+^2 ^^3 „ 

e 2 +e 2 > e is a suffieient one. 


The asymptotie analysis ean be performed as whose details are omitted here. How¬ 

ever, it should be remarked here that two-soliton solution for all eomponents in this ease always un- 
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dertakes elastic collision without shape changing. This feature is the same as the one of the mixed 
soliton solution including one-bright-two-dark soliton for SW components in 2D ( 3 - 1-1 )-component 
YO system.^^ In Fig. 6 , we exhibit this mixed-type soliton solution under the same nonlinearity co¬ 
efficients with Fig. 5 and the parameters are given as pi = P2 = OCi = 1, OC2 = 2 ,/>i — —^—i,p2 — 
^ -I- i, = 1 -I- i, = 2 -I- i and 0io = 020 = 0- Fig- 6 (a) shows two bright solitons collisions for 
the SW components 5 *^^^ pig 5 (p) and (c) display two dark solitons collisions for the SW com¬ 
ponents and with the different amplitudes, and Fig. 6 (d) represents two bright solitons 
collision for the LW component —L. 

For the construction of the bound states from the mixed soliton solution (59)-(61), one can 
derive the same condition pu = Pii as the previous subsection. By taking the same nonlinearity 
coefficients in Fig. 5, such a bound state with the parameters Pi = p 2 = OCi = 1 , 0 C 2 = = 

— \ — = ^ — = 1 ,C 2 ^^ = ^ and 0io = 02 O = 0- is depicted in Fig. 7. Fig. 7 (a) exhibits 

two-bright-soliton bound state for the SW components ^ pig. 7 (P) and (c) display two-dark- 
soliton bound states for the SW component and with different amplitudes and two-bright- 
soliton bound state of the LW component —L is shown in Fig. 7 (d). 

Here we need to point out that A^-soliton bound state can be formed in both above mixed 
soliton solution for arbitrary combination of the nonlinearity coefficient (£ = 1,2,3). More 
specifically, if the wave number for the /-th soliton among an N-soliton solution is pt = pm -f ipn 
, then Qm = piR{x + Iput) - 1 - 0 ,which suggests that the velocity for the i-th soliton is —Ipu. 
Thus a multiple bound state in which solitons move with the same velocity only requires all the 
Pu to be the same. However, higher-order soliton bound state cannot exist in the soliton solution 
comprised of all dark soliton for the SW components. As reported in,^^ due to the existence of the 
parameter constraint in general multi-dark soliton solution, only two-dark-soliton bound state is 
possible under the condition that the nonlinearity coefficients take different signs. 



Fig. 6. (Color online) Mixed two-soliton solution (one-bright-two-dark soliton for SW components) in the (3+1)- 
component YO system. 
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Fig. 7. (Color online) Mixed two-soliton bound state (one-bright-two-dark soliton for S W components) in the (3 h- 1)- 
component YO system. 


3. General mixed soliton solution to the ID multicomponent YO system 

In the same spirit as the (3-1-1 )-component YO system, the general mixed type soliton solution 
to the ID {M + l)-component YO system can be derived by the KP hierarchy reduction method. 
It is known that the multi-bright soliton solutions can be derived from the reduction of the multi- 
component KP hierarchy, whereas, the multi-dark soliton solutions are obtained from the reduction 
of the single KP hierarchy but with multiple copies of shifted singular points. Therefore, if we con¬ 
sider a general mixed soliton solution consisting of m bright solitons and M — m dark solitons to the 
SW components in the ID (M-f l)-component YO system (1), we need to start from an (m-f 1)- 
component KP hierarchy with M — m copies of shifted singular points in the first component. 
Furthermore, by performing the complex conjugation and dimension reductions, these bilinear 
equations become ones of multicomponent YO system. Meanwhile, the general bright soliton so¬ 
lution can be reduced from the tau functions of the KP hierarchy. The details are omitted and we 
only provide the final result here. 

To seek for mixed multi-soliton solution consisting of m bright solitons and M — m dark ones 
for SW components, the ID multicomponent YO system is first converted to the following bilinear 
form 

[iA-Dgg(^)-/ = 0, k=l,2,---,m, 

[i(A-2a/D;,)-DgfiW-/ = 0, I = 1,2,-■ ■ ,M-m, (63) 

(AD, - 2 0/+„p2) / ■ / = -2Er - 2 LZT , 

through the dependent variable transformations: 

= = pz^e‘(“'^+“/0, L=-2(log/),„ (64) 

where a/, |3/ and p/(p/ > 0) are arbitrary constants. 

Similar to the procedure discussed in Section 2, one can obtain mixed multi-soliton solution as 
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follows: 


/ = 


A 1 
-I B 


, 


A« 

-7 


1 

B 




A 7 4)^ 
-7 B 0^ 


0 a 


where A, and 5 are N xN matrices whose entries are 


1 


Clij — 


Pi + P*i 


-e®'+®ri 


(0 _ 1 

A + 7?' 


bij = 




Pi - IOC/ 

P* + ioc/ 


iEJLiOtcf’Vf 


(p* + Pj) b,’ -Pj+z1i" (,,.+i„,)(pj-i„,)j 
and (|) and Q are /7-component row vectors given by 


i^/+mP/ 


♦ =(e«'.e''S...,e"»). Q =-(cl-.c™.. ■. ,c™). 

with 0, = — ipjt + 0,0 and pi, 0,o and cf’ {i = 1,2, ■ ■ ■ ,N) are arbitrary complex parameters. 

In the above solution, one necessary condition similar to the (3-1-1 )-component YO system for 
the existence of an /7-soliton solution is found as follows 


(k) (k) 


Sk)-s 


0 


(65) 


E C)fc|cJ 


ik )\2 


M—m 




( 66 ) 


As reported in,^^ the arbitrariness of nonlinearity coefficients 0 £ increases the freedom resulting 
in rich mixed soliton dynamics. Here the present solution admits mixed A-soliton in the ID mul¬ 
ticomponent YO system for all types of nonlinearity coefficients, including positive, negative and 
mixed types. 

The formula for general mixed soliton solution can be generalized to include all bright and 
all dark soliton solutions, which is the same as the vector NLS equation.^^ More specifically, the 
general bright soliton solution can be viewed directly as a special case of the bright-dark soliton 
one. When m = M, the A-bright soliton solution takes the same determinant form as the general 
bright-dark one. The general dark soliton solution was also derived from the single KP hierarchy 
but the corresponding dimension reduction results in the distinct parameters constraint. The final 
expression of the A dark-soliton solution is different from the one of the bright-dark soliton solu¬ 
tion. However, when m = 0, as discussed in [38], an alternative form of all dark soliton solution 
takes the same form as the solution (65) by redefining the matrix B as an identity matrix (Bij = dij) 
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and imposing the parameters constraint as follows 


p*i-Pi+Y. 


(67) 


;“i 

Due to a simple determinant identity shown in,^^ this form of dark soliton solution coincides with 
the one in.^^ 


4. Summary and conclusion 

We have constructed the general bright-dark A^-soliton solution to one-dimensional multicom¬ 
ponent YO system describing the nonlinear resonant interaction of M short waves with a long 
wave, i.e., one-dimensional (M + l)-component YO system. This solution exists in the original 
system for all possible combinations of nonlinearity coefficients including positive, negative and 
mixed types. Taking the (3-1-1 )-component YO system as an example, we have deduced two kinds 
of the general mixed A^-soliton solution (two-bright-one-dark soliton and one-bright-two-dark one 
for SW components) in the form of Gram determinant by using the KP hierarchy reduction method. 
Then, the same analysis was extended to obtain the general mixed solution consisting of m bright 
solitons and M — m dark ones for SW components in the {M + l)-component YO system. The 
expression of the mixed solution also contains the general bright and dark A^-soliton solution. 

For the dynamics of the mixed solitons, in parallel with the 2D multicomponent YO system,^^ 
the energy exchanging collision of mixed solitons can be realized only in the bright parts of the 
mixed solitons and is possible only if the bright parts of the mixed solitons appear at least in 
two SW components. Such an interesting phenomenon can also be found in ID (3-1-1 )-component 
YO system as discussed in previous section, as well as in the ID multi-component YO system. 
In addition, the related analysis regarding mixed-soliton bound state in ID (3-1-1 )-component YO 
system implies that arbitrary higher-order mixed-soliton bound state in ID multi-component YO 
system can also be formed. 
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